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Abstract 

We report experimental observation of the conversion of a phase-invariant nonlinear system into 
a phase-locked one via the mechanism of rocking [de Valcarcel and Staliunas, Phys. Rev. E 67, 
026604 (2003)]. This conversion results in that vortices of the phase-invariant system are being 
replaced by phase patterns such as domain walls. The experiment is carried out on a photorefractive 
oscillator in a two-wave mixing configuration. A model for the experimental device is given that 
reproduces the observed behavior. 
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Introduction— The symmetry properties of the order parameter phase in extended non- 
linear systems determine the type of localized structures supported by these systems. On 
the one hand phase-invariant systems display vortices, which are phase defects of the or- 
der parameter, around which the phase changes in 2tt, forcing the order parameter to be 
null at the vortex center. Self-oscillatory systems represent a paradigm of such behavior 
as the oscillating state reached after a homogeneous Hopf bifurcation can have any phase 
(the system is autonomous). Examples of such systems are certain chemical reaction (like 
the Belousov-Zhabotinsky reaction) and several optical systems (like lasers, nondegenerate 
optical parametric oscillators, OPOs, and nondegenerate wave mixing photorefractive oscil- 
lators). On the other hand there are nonlinear systems with broken phase invariance that, 
obviously, cannot support vortices. Among them, especially interesting systems are those 
displaying phase bistability. Such systems support a different type of localized structure, 
the domain wall, which connects spatial regions where the orderparameter passes from one 
homogeneous state, F , to the equivalent, symmetric state, — F |l| so that the phase changes 
by 7r from one side of the wall to the other one. A paradigm of such type of phase-bistable 
system is the degenerate OPO. 

The question we address here is: Is it possible to convert a phase-invariant system into a 
phase-bistable one by some simple external action? In other words, can a system displaying 
vortices be forced to display domain walls? An old and well known answer to these questions 
consists in the periodic forcing of a self-oscillatory system at a frequency around two times its 
natural oscillation frequency uj$ (resonance 2:1). Under that type of (so called) parametric 
driving the originally phase-invariant system is predicted "\ and experimentally observed 0, 
Q to transform into a phase-bistable system, exhibiting domain walls. But while parametric 
driving is useful in many contexts it is not so in general in nonlinear optics. For instance, 
a laser emitting at frequency Uq is insensitive to a forcing at 2u; because the gain line is 
extremely narrow as compared with the magnitude of the optical frequency ujq. Yet, a new 
mechanism coined rocking [5] has been recently proposed to overcome such limitation of the 
usual parametric driving. Rocking consists in forcing a self-oscillatory system around its 
natural oscillation frequency (resonance 1 : 1) so that the forcing amplitude is periodically 
modulated in time at a frequency uo <C uu . Thus rocking is a multifrequency forcing (in its 
simplest case, a sinusoidal modulation of the forcing amplitude, it is a bichromatic forcing) 
around the resonance 1 : 1 of a self-oscillatory system. We want to remark that unlike the 
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parametric driving, any self-oscillatory system (including lasers and any nonlinear optical 
system) should be sensitive to rocking as this forcing acts on the main resonance of the 
system. 

In order to gain insight into the rocking idea, consider a particle in a 2D space with 
coordinates (x, y) under the action of a potential with the sape of a sombrero -having a 
maximum at the origin (0, 0) and a degenerated minimum at x 2 + y 2 = r^ in . This is a phase- 
invariant system, what means that the equilibrium position is phase(angle)-degenerated. 
Suppose now that this potential is rocked along the x-direction (e.g., by adding a term 
x sin ut to the original potential). In this case the phase degeneracy is obviously broken 
and now the system tends to be around (x = 0, y — ±r min ): The system is now phase- 
bistable. This pictorial image illustrates the physical rationale behind the rocking idea: 
The coordinates of the fictious particle correspond to the real and imaginary parts of the 
complex amplitude of oscillations of the self-oscillatory system, whose evolution derives from 
a sombrero-like potential in the simplest case, as originally introduced in 5} through an 
analysis of a complex Ginzburg-Landau equation, which is the simplest model for spatially 
extended self-oscillatory systems [jj. As shown in j^j, under the action of rocking the order 
parameter of the system F(r,t) (e.g., the laser electric field complex amplitude) develops 
two components, F (r, t) = F^ (t) (r, t), where F^ is a ^j-periodic function of time that 
follows the modulation, at frequency u, of the forcing amplitude and ip is a phase-bistable, 
reduced order parameter that can display domain walls, phase domains, and phase domain 
solitons 5]. Even if the original prediction is based on a complex Ginzburg-Landau model, 
rocking is expected to be of wide applicability as the phase-bistability mechanism introduced 
above seems to be quite model independent. 

In this Letter we give the first experimental evidence of rocking induced phase-bistability 
and the associated formation of phase domains and domain walls. Our system is a (nonde- 
generate) two- wave mixing photorefractive oscillator (PRO), which highly resembles a laser 
from the nonlinear dynamics viewpoint. Rocking is done by injecting an amplitude modu- 
lated laser beam into the resonator. We show that in the absence of rocking (free running 
configuration) the PRO exhibits vortices, which are converted into phase domains under 
the action of rocking. Experiments performed under a quasi ID transverse geometry (in 
order to avoid curvature effects) show that the system displays domain walls in a window 
of rocking modulation frequencies, in agreement with 5]. A model for the experimental de- 
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vice is presented that reproduces the experimental findings. In a limit, such model reduces 
to a parametrically driven complex Swift-Hohenberg equation, which represents a general- 
ization of to other types of order parameter equations and hence enlarges the range of 
applicability of rocking. 

Experimental setup - The PRO, Fig. 1, consists of a photorefractive BaTi03 crystal 
placed inside a (near) self-imaging ring resonator [6| of cavity length 1.2m (the resonator 
free-spectral range is 250MHz), similar to that used in [7(. The effective cavity length is 
approximately 2cm, which is actively stabilized by means of piezomirror PZT1 in order 
to have a precise control on cavity detuning (the difference between the frequencies of the 
pump field and that of the closest cavity mode). The crystal is pumped by a singlemode 
514nm Ar + laser with a power around lOOmW cm -2 . An amplitude modulated beam (the 
rocking beam) coherent to the pump field is injected into the cavity. The c+ crystal axis is 
oriented so that gain is maximized when all three beams (pump, oscillation and injection) 



are extraordinarily polarized. The intracavity slit D is p 
make the system quasi ID in the transverse dimension 



aced in a Fourier plane in order to 
Finally, there is the possibility 



of injecting a tilted coherent beam in order to "write" domain walls as in 0, lflj ]. 

The rocking beam must be an amplitude modulated field of zero mean: We chose to 
use a field of constant amplitude, whose phase changes exactly by 7r in every half period. 
We do this by injecting into the cavity a beam coming from the pumping laser, after being 
reflected on piezomirror PZT2, Fig. 1, which is moved periodically back and forth by half a 
wavelength (the modulation frequency is on the order of 1Hz). This way the rocking has a 
pure amplitude modulation (only the sign of the field amplitude changes). This operation 
must be very precise as phase jumps different to n do not produce the desired result. 

Experimental results. - The cavity length is chosen for the system be in (almost) exact 
resonance. This is not mandatory but is the simplest way for the rocking field be resonant 
with the intracavity field (see below the model section). When the rocking beam is off and 
the intracavity slit is open, the system spontaneously forms vortices in the output field, as 
expected U|. Under these conditions, the application of a rocking beam is able to transform 
vortices into phase domains, as shown in Fig. 2. This result is a direct demonstration of the 
rocking induced phase-bistability. 

Two dimensional phase domains are transient structures due to curvature [ij. Hence we 
performed a series of experiments under quasi ID conditions (by narrowing the intracavity 
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slit), which allow stable domain walls in phase bistable systems |9[. These domain walls are 
however unstable in ID phase invariant systems: If a domain wall is injected into the free 
running PRO, Fig. 3(a), we observe that the domain wall vanishes with time, Fig. 3(b), and 
is replaced by a spatially uniform state, Fig. 3(c). On the contrary, if the PRO is rocked 
completely different results are obtained. Now an injected domain wall, Fig. 3(e), remains 
stable and fixed, Figs. 3(f) and (g). (The amplitude modulation frequency of the rocking 
beam was 1.5Hz and its intensity was of the same order as that of the output field in the 
absence of rocking.) Figure 3(h) evidences that we are in the presence of an Ising wall 
The structure is robust and this proofs that the system is now phase bistable. The ability 
of rocking to sustain domain walls is thus proven. 

There remains to clarify which is the effect on the performance of the system of changing 
the modulation frequency and amplitude of the rocking beam. We have observed that the 
frequency of modulation cannot be too small or too large: For modulation frequencies of 
0.1Hz or below rocking cannot sustain the injected domain wall. The same happens for 
modulation frequencies greater than 10Hz. It seems that modulations from 1 to 3Hz are 
optimal, which are on the order of the inverse of the photorefractive grating decay time. 
While we do not have definite measurements, we can say that rocking is effective in a 
window of rocking intensities as well. 

All previous featu.es a r e to a„t with the pactions of Re, g The theory is 
lowever based on a Ginzburg- Landau model, which should not be valid near resonance 
the ease eonsidered here. In order to justify theory the r esu.ts ou a m0r e 
firm basis we thus proceed to model the PRO. 

Theory.- We adopt the two-wave mixing PRO model for a purely diffusive photorefractive 
crystal |l2j, such as BaTiOs, generalized to account for the injection of the rocking beam 
(details will be given elsewhere). The model equations, suitably normalized, can be written 
as: 



a-%F = - (1 + iA) F + iaV 2 F + N + F R {t) , (1) 
d t N = -N + g— F —^, (2) 

where F (r, t) is the slowly varying envelope of the intracavity field, N (r, t) is the complex 
amplitude of the photorefractive nonlinear grating, r = (x, y) are the transverse coordinates, 
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V 2 = d 2 + dy, a = kt is the product of the cavity linewidth k with the photorefractive 
response time r (a > 10 8 under typical conditions, and r ~ Is), t is time measured in units 
of r, the detuning A = (u c — up) / n (up and u c are the frequencies of the pump and of its 
nearest cavity longitudinal mode, respectively), a is the diffraction coefficient (that depends 
upon geometry [8( and can take either sign), Fp is the complex envelope of the rocking 
(injected) field, and g is the (real) gain parameter that depends on crystal parameters and 
on the geometry of the interaction. The actual intracavity field, S, and rocking field, Sp, 
read £ = Re FEpe~ luJprt and Sp = Re FpEpe~ luJpTt , respectively, where Ep and up are the 
complex amplitude and the angular frequency of the pumping laser field. Note that Ep acts 



just as a scaling factor (it is unrelated with the gain parameter g) 121. Il3|. 



For Fp = the model is equivalent to that in [12j , which holds the continuous symmetry 
(F, N) — > (Fe 1 ^, Ne^) . Hence the system is phase invariant in the absence of forcing. 
This free running PRO model has two main solutions jl2j ]: The trivial solution F = N = 
0, and the family of traveling- wave solutions (parametrized by the wavevector k) F = 
^/g/g - le^-'-W and N = (1 + ifl) F, with g = and Q = ^ (A + ak 2 ) -> A+ak 2 

as a 3> 1. One easily sees that: (i) for A/a > the oscillation threshold is minimum for 
k = (on-axis emission), occurs at g = 1 + A 2 , and the frequency of the generated field 
is shifted by Q = A from that of the pump beam; and (ii) for A/a < the threshold is 
minimum for k = a/— A/a (off-axis emission), occurs at g — 1, and there is no frequency 
shift (fl — 0). We checked these features in our experiment. In particular cavity resonance 
was determined by interfering the emitted field with a reference coming from the pumping 
laser: The cavity length at which the beating frequency JQ) passes from zero to a non 
null value (or vice versa) corresponds to exact resonance |l2j]. These facts, together with 
the phase invariance of the model, make the PRO in the two-wave mixing configuration a 



system largely equivalent to a large aspect ratio laser |12|, |14 1 . 

In the following we use Fp = R cos ujt (R real without loss of generality), which is the 
simplest form of amplitude modulation and corresponds to a bichromatic injected signal. In 
order to give analytical evidence of the rocking induced phase-bistability in the PRO model 
we perform an asymptotic expansion of Eqs. (jll2J) based on the of multiple scales technique 
[15I ] . In order to approach the experimental conditions we assume small detuning, and large 
cavity linewidth: We take A = O (e) , a = O (£ -4 ) with e a smallness parameter (0 < e 1) 
(the final result does not depend on the precise scaling for cr, whenever o ^> 1 is assumed). 
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We further assume that gain is close to threshold, g = 1 + 0(e 2 ), and that the rocking 
jarameters verify R — O (e), u — O (1). The analysis, similar to that performed in laser 



16| and nondegenerate 0P0 |l7|, [l8| models, yields N = [1 + i (A — aV 2 )] F — Rcosout, 



F(r,t) = F w (t) + itlj(r,t) + 0(e 3 ) J (3) 



the order parameter r/> verifying a 



jr,(t) = J (1 - 2iA)sinu;t + (^1 - zA^^-j cjcoscjt 
complex Swift-Hohenberg equation with parametric gain: 

<9 t v = 7</>* + G? - i - 27) <A - M 2 

+ i (aV 2 — A) -0 — (aV 2 - A) 2 ?/>, (4) 

where 7 = "y-^r- is the rocking parameter (note that 7 > ; y-). In the absence of rocking 
(7 = 0), Eq. (j3J) is phase invariant and is isomorphic to those describing lasers lid and 
nondegenerate OPOs |l3,I3], as well as highly resembles that for drift-type PROs The 
role of rocking is clearly appreciated in Eq. (j3J): It introduces an phase-sensitive gain (first 
term in the r.h.s.) that breaks the original phase invariance of the undriven system down 
to the discrete one ip — > —ip. Thus the PRO becomes phase bistable. The spatially uniform 
steady states of Eq. (J1J) are given by ip = ± e llfi , \ip\ 2 = fi — 27 + A/7 2 — A 2 , \i = g—1 — A 2 , 
e 2ttp = ^ 1 ~^ — (other two, intrinsically unstable states exist as well, which we do not 
consider). These phase-locked states exist if |A| < 7 < 7 max (7max = ^ + |\//i 2 — 3A 2 ) as 
far as /i > v^3 |A|. These inequalities imply analogous ones for the rocking intensity, i? 2 , or 
the rocking modulation frequency, u, depending on which parameter (u or R) is kept fixed. 
This prediction supports the experimental findings described above. Finally note that as 
7 > -y- the rocked states will exist only if R 2 < 27 max , and [i > \/3|A|, for any to. 

In a series of numerical experiments we have checked the above predictions and have 
found that they remain valid even far from the asymptotic limit described by Eq. (@J), as 
it happens in the original proposal jsj]. In order to have additional comparison between 
theory and experiment, we present in Fig. 4 some numerical simulations of Eqs. (|1I2() for 
a = 10 2 (large, but not extremely large in order to avoid stiffness problems), A = 0, g = 2 
(we estimate that the gain in the experiment is about 100% above threshold), u = 2tt and 
R = 0.5. Comparison with Fig. 3 shows that results are very similar to the experimental 
ones. A message from this theoretical treatment arises: Phase invariant systems described 
by order parameter equations of different nature (like the Ginzburg-Landau and the Swift- 
Hohenberg complex models) behave similarly under the influence of rocking. 



In conclusion, we have experimentally demonstrated, and theoretically justified, the abil- 
ity of the rocking mechanism introduced in U to generate phase-bistable states in otherwise 
phase invariant systems. The studied system, a PRO in two-wave mixing configuration, 
highly resembles laser systems from the nonlinear dynamics viewpoint. Hence the results 
put forward in this Letter should motivate similar experiments in laser systems, which could 
have potential applications in the field of information technologies. 
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discussions with Javier Garcia Monreal and Kestutis Staliunas. 
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Figure captions 

Fig. 1.- Scheme of the experimental setup. A photorefractive BaTi03 cristal is placed 
inside the quasi self-imaging resonator formed by mirrors Ml, M2 and M3, piezomirror PZT1 
and the four lenses of focal length /. D is a diafragm, and PZT2 is a piezomirror driven by 
a square-wave signal, plot V vs. time, that produces the rocking beam. 

Fig. 2.- Interferometric snapshots of vortices (a) existing in the free running PRO and 
phase domain structure (b) generated in the rocked PRO at resonance. Magnifications of 
the interferometric images in the position of (c) a vortex showing the annihilation of the 
interference fringes and (d) of a part of the domain wall showing the ir phase jump. 

Fig. 3.- Experimental snapshots of injected domain walls (DW) close to cavity reso- 
nance. (a)-(c): Disappearance of a DW in the free running PRO (the time interval between 
snapshots is 5s). (e)-(g): Stabilization of the DW in the rocked PRO (the rocking frequency 
is 1.5Hz and the time interval between snapshots is 15s). (d) and (h): Horizontal cuts show- 
ing the field amplitude (a.u.) and phase corresponding to snapshots (b) and (f), respectively. 
The snapshots transverse dimension is 1.2mm. 

Fig. 4.- As Fig. 3 but from numerical simulations of Eqs. (0), (J2J). Parameters: a = 10 2 , 
A = 0, g = 2, uj = 2tt, and R = 0.5. The length of the horizontal dimension is 21-y/a in 
(a)-(c) and (e)-(g), lk\/a in (d) and 9.3y/a in (h). 
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Figure 1 

A. Esteban-Martin et. al, Rocking of a nonlinear optical. 



Figure 2 

A. Esteban-Martin et. al, Rocking of a nonlinear optical... 




Figure 3 

A. Esteban-Martin et. al, Rocking of a nonlinear optical. 
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Figure 4 

A. Esteban-Martin et. al, Rocking of a nonlinear optical. 



